Charge-order transition in the extended Hubbard model on a two-leg ladder 
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We investigate the charge-order transition at zero temperature in a two-leg Hubbard ladder with 
additional nearest-neighbor Coulomb repulsion V using the Density Matrix Renormalization Group 
technique. We consider electron densities between quarter and half filling. For quarter filling and 
U — St, we find evidence for a continuous phase transition between a homogeneous state at small 
V and a broken-symmetry state with "checkerboard" [wavevector Q = (n, n)] charge order at large 
V . This transition to a checkerboard charge- ordered state remains present at all larger fillings, but 
becomes discontinuous at sufficiently large filling. We discuss the influence of U/t on the transition 
and estimate the position of the tricritical points. 



The competition between kinetic energy and Coulomb 
repulsion in electronic systems can lead to a variety of 
interesting phenomena, one of them being charge order- 
ing. A periodic charge order, i.e., a modulation of the 
charge density, can be described as charge density wave. 
One possible mechanism for charge ordering is the crys- 
tallization of electrons due to their long-range Coulomb 
repulsion as proposed by Wigner.El The Wigner lattice of 
electrons forms without the underlying lattice structure 
at low densities when the system is dominated by the 
effect of the Coulomb repulsion. Charge ordering may 
also occur at higher electron densities if the kinetic en- 
ergy is reduced due to small hybridization of orbitals, 
or due to the interaction with lattice or spin degrees of 
freedom. Charge-ordered states have been observed in, 
for example, rare earth manganites, which have attracted 
attention recently due to their "colossal" magnetoresis- 
tance. Several of these compounds (e.g. Lai_a;Ca K Mn03 
for x > 0.5) show a charge-ordered ground state for a cer- 
tain range of doping.B Another material showing charge 
order is NaV2 05. It undergoes a phase transition at 
T c = 34 K that is characterized by the opening of a spin 
gap and a doubling of the unit cell. While this tran- 
sition was originally thought to be spin-PeieeJs, recent 
studies have found evidence for charge order. 013 It has 
been proposed that this material is well-described as a 
quarter-filled ladderBEl Two-leg ladder models are also 
thought to be relevant to a number of other materials 
containing ladder-like structures, such as the the Vana- 
dates MgV205 and CaV20s, and the cuprates SrCu203 
and Sri4Cu2404i. For a more detailed description of lad- 
der materials and models, we direct the reader to a recent 
reviewa and the references contained therein. 

One of the simplest models of interacting electrons that 
allows for charge ordering is the Hubbard model extended 
with an additional nearest-neighbor (NN) Coulomb re- 
pulsion, V. The charge order transition in this model 
has been studied in the one-dimensional (ID) model in 



dimensions) at quarterly and half filling.Ej A variety of 
techniques, such as mean-field approximations, perturba- 
tion theory, and numerical methods as Quantum Monte 
Carlo and the Density Matrix Renormalization Group 
(DMRG) have been employed in these studies. Their re- 
sults can be summarized as follows: At the mean-field 
level, the transition between a homogeneous state and a 
charge-density wave (CDW) state at half filling in a hy- 
percubic lattice occurs at V c = U/zq where zq denotes 
the nupber of nearest neighbors (zq = 2d). Numerical 
studiesE3 indicate a slightly higher value of V c , at least 
in ID. Interestingly, in ID at half filling the transition 
has been found to be second order at small U/t and first 
order at \&rgeJJ/t with the tricritical point located at 
U c /t ~ 4 — 6.E£l Here we use the term "first order" to 
denote discontinuous behavior of the charge order pa- 
rameter as a function of microscopic parameters such as 
V or band filling, and "second order" to denote contin- 
uous behavior. At quarter filling in ID, the situation is 
more complicated since a number of phases compete for 
large VS3 In higher diajiaisions, the transition seems to 
always be first order.EJila However, conclusive studies 
that can reliably distinguish between first- and second- 
order transitions are lacking. Between quarter and half 
filling, the existence and nature of a transition is unclear 
in general. In ID, the extended Hubbard model under- 
goes phase separation rather than_a transition to a CDW 
state for \U\/t < 4 and large Vt3 For U/t > 4, indica- 
tions are that the ID system undergoes a transition to_a 
q = n CDW state for sufficiently large V at all fillings.E£l 
In this paper, we examine the charge-order transition 
in the extended Hubbard model on the two-leg ladder, 
considering all values of band filling between quarter and 
half filling. We shall investigate the nature and location 
of the charge-order transition in the ladder and compare 
with the ID as well as higher dimensional models. 

The single-band extended Hubbard model has the 
Hamiltonian 
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Here we consider a lattice consisting of two chains of 
length L, i.e., a ladder, and discuss band fillings (n) — 
N/(2L), with TV the number of electrons. The summa- 
tion (ij) then runs over all pairs of nearest neighbor sites 
in the ladder. In this work, we take both the hopping and 
the nearest-neighbor Coulomb repulsion to be isotropic; 
in general, one could consider the anisotropic case, pa- 
rameterized by tii, t±, Vji and V±. 

The numerical results shown in this-.paper have been 
calculated with the DMRG techniquecJ on lattices of up 
to 2 x 64 sites with open boundary conditions both be- 
tween the two chains and at the ends of the chains. We 
have kept up to 600 states per block, resulting in the 
discarded weight of the density matrix eigenvalues being 
typically 10~ 8 or less. The errors in the energies and 
correlation functions arising from the truncation of the 
density matrix are always less than a few percent. 



which is defined similarly in terms of the spin-spin cor- 
relation function (Sj +i Sj). 

In Fig. [j], we show the results for C(q) and 5(q) for 
a 2 x 64 system at U = 8t and quarter filling for dif- 
ferent values of V. The transition between a homoge- 
neous state at small V and an ordered CDW state at 
large V can clearly be seen. The DMRG ground state at 
V > V c is fully gapped and has broken translation invari- 
ance, i.e., for large V every second site is empty. [This 
gapped charge-ordered state with (jr, it) broken symme- 
try is present at all fillings from quarter filling to half 
filling.] In the spin channel, peaks at (0, w) and (w, 0) de- 
velop with increasing V showing an alternating spin pat- 
tern on the occupied sites in the CDW state. This spin 
order arises from virtual hopping processes (fourth order 
in t) which lead to antiferromagnetic couplings between 
the occupied sites. In the limit of large V the system 
can be mapped to a J\ — Ji spin- \ chain where J\ ^> Ji 
because of the larger number of exchange processes in 
diagonal (1,1) direction. 
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FIG. 1. Charge- and spin correlation functions for a 64 x 2 
ladder, U = 8i, quarter filling (n) = 0.5 (64 electrons), and 
different values of the nearest-neighbor repulsion V. 

We now turn to the discussion of the results. In order 
to investigate the charge ordering, we consider the static 
charge structure factor 



where 



C(q) = i^e^ Ri C(R i ) 



av r -i 

in 



(2) 



(3) 



(...) denotes the ground-state expectation value, Srij — 
rij — (rij), and we average over typically N av = 6 sites 
to remove oscillations due to the open boundaries. For 
checkerboard charge order, which is expected to be the 
ground state for large values of V, C(q) should show 
a pronounced peak at Q = (ir,ir). The spin order can 
be studied by looking at the spin structure factor, S(q), 
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FIG. 2. Finite-size scaling for the staggered charge cor- 
relation function at U = St and (n) = 0.75. The curves cor- 
respond to V/t = 2.5,2.8,2.9,2.95,3.0,3.05,3.1,3.2,3.5,4.0, 
and 8.0 from bottom to top. The dashed lines are linear fits 
through the data points for L = 32, 48, 64. The irregular 
behavior at V/t — 3.05 may indicate a first order transition 
(here the L — 32 point was excluded from the fit). 

At half filling, the spin correlations correspond to 
short-range (it, ir) order at small V/U. In the limit of 
large V 3> (t, U) the spin correlations vanish. This be- 
havior is due to the increase in the number of double 
occupancies in the charge-ordered state at half filling, 
which leads to a suppression of local magnetic moments. 
Between quarter and half filling, 5(q) shows peaks at in- 
commensurate wavevectors that correspond to combina- 
tions of bonding and antibonding Fermi vectorso. Near 
quarter filling, these peaks tend to be strengthened but 
also displaced as V is increased, whereas closer to half fill- 
ing the peaks tend to be suppressed in the CDW phase. 
We further note that between quarter and half filling the 
residual magnetic moments tend to order ferromagncti- 
cally in the charge-ordered phase, i.e., the ground-state 
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spin is non-zero.u The existence-pf gapless modes in the 
homogeneous phase is less cleared; at least at half filling 
both spin and charge gaps are non-zero for any parame- 
ters with U > 0. 

In order to examine the extent of the charge ordering 
and the nature of the transition, we calculate the order 
parameter for a Q = (n, n) charge-ordered state, 
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We carry out the L — > oo extrapolation by calculating 
r)(L) for L — 16,32,48,64 and performing a linear fit 
in 1/L through the three largest system sizes, see Fig. 
0. (The L = 16 point is consistent with the linear ex- 
trapolation except in the region of the transition.) This 
extrapolation is the major source of uncertainty in our 
results; the error here is typically A77 = 0.01. 

In Fig. g, we plot y/rj (which corresponds to the 
relative difference of the sublattice occupancies in the 
broken-symmetry charge-ordered state) for U = 8t and 
different band fillings. For U = 4£ as well as U — 8t, we 
find a transition from a homogeneous state to a charge- 
ordered state with increasing V for all values of (n) . Note 
that the extrapolation with system size is crucial in or- 
der to extract the V c from our data. Examination of 
rj(V;L) for fixed L, as used in Ref. [12| to determine V c 
from DMRG results for the ID system, is not adequate 
for the ladder system. 




FIG. 3. Order parameter ^/rj versus V/t for U = 8t and 
(n) = 0.5, 0.75 and 1. The data indicate a second-order 
quantum phase transition at V c /t = 2.65 ±0.1 between a 
homogeneous and a charge-ordered state for (n) = 0.5. For 
(n) = 0.75 the transition at V c /t — 3.05 ±0.1 is weakly first 
order whereas the transition at half filling at V c /t = 2.8 ±0.1 
is clearly first order. 

An additional, related method to determine V c is to 
examine the scaling of L ■ C(Q; L) with L. This quan- 
tity should diverge with L in the charge-ordered state, 
but should decrease with L in the homogeneous state to 
a value determined by the density correlation length of 
the system.Eil This method yields the same values of V c 
as that obtained from the extrapolation of y/rj to within 



A V c « O.li. Furthermore, we have checked that the com- 
pressibility of the system is positive, i.e., there is no ten- 
dency to phase separation in the parameter region stud- 
ied here (U/t > 4). Note, however, that phase separation 
occurs for smaller values of U/t, e.g., for U/t = 1 and 
V ^> t. A complete investigation of the phase diagram 
will be published elsewhereO 

At U = 8t (Fig. ||) the transition to the charge-ordered 
state is continuous for small filling (e.g. quarter filling, 
(n) = 0.5) whereas for large filling the order parameter 
clearly shows a jump at a critical value V c . This leads 
to a tricritical point in the V—(n) plane, as can be seen 
in the phase diagram of Fig. 0. This tricritical behav- 
ior seems to be similar.-tp the one observed in the ID 
system at half fillingO'Ej The location of the tricritical 
point depends on the strength of U, as shown in Fig. |j. 
For U = 4£, its position is (n) tc = 0.95 ± 0.1, and at 
U = 8t we find (n) tc = 0.65 ± 0.1. The position of the 
tricritical point (for fixed U) has been estimated from the 
(n)-dependence of the magnitude of the jump in i](V). 
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FIG. 4. Phase diagram in the V-(n) plane showing first 
order (solid lines) and second order (dashed lines) transitions 
between the homogeneous and the charge-ordered state for 
U = At and U — 8t. The boundaries have error AV C ~ O.lf 
which arises mainly from the L — > 00 extrapolation. The cir- 
cles indicate the estimated positions of the tricritical point 
with uncertainty A{n tc ) ~ 0.1. 

Finally, we address the dependence of V c on U and 
the band filling. At halfJiUing, both weak- and strong- 
coupling approximations ! 1 6 tP 2 l for the transition between 
spin-density wave (SDW) and CDW yield a critical value 
of V c — U I zq for a system on a hypercubic lattice. How- 
ever, these Hartree-type approximations assume a SDW 
phase with true long-range order which is absent for ID 
chains as well as for the ladder system considered here. 
Nevertheless, the mean-field V c is in reasonable agree- 
ment with the value of slightly more than (7/2 found nu- 
merically for a ID chainJij For the ladder systems studied 
in the present work, the DMRG results for half filling also 
agree quite well with the predictions of the Hartree ap- 
proximation. As can be seen in Fig. || U/V c is consistent 
with V c = (7/3 in both the weak- and strong-coupling 
limits. 



Below half filling, the behavior is more complicated 
because it is no longer dominated only by the interplay 
of U and V. Since the average number of double oc- 
cupancies decreases with decreasing electron density, the 
kinetic energy becomes important for larger U. In the re- 
gion between quarter filling and half filling, we are faced 
with the interplay of all three energy scales t, U, and V, 
and as can be seen in Fig. |J, V c varies considerably with 
(n) and U. 
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FIG. 5. The ratio U/V c at half filling for different values of 
U. The data indicate that U/V c approaches the value zo = 3 
for both weak and strong coupling which is consistent with 
mean-field arguments. 

In summary, our scenario for the charge-order transi- 
tion in the Hubbard ladder with U > At proceeds as fol- 
lows: At all band fillings, an increasing nearest-neighbor 
repulsion V drives a transition from a homogeneous state 
to a gapped CDW state with Q = (tt, it). At quarter fill- 
ing, each occupied site in the CDW state carries spin |, 
and the residual kinetic energy orders these spins "anti- 
ferromagnetically" at ordering vectors (n,0) and (0,7r). 
A similar charge ordering leading to "zigzag" antiferro- 
magnetic ordering has been found at quarter filliag in 
coupled Hubbard ladders within mean-field theoryQ and. 
in the t-J-V ladder with adiabatic Holstein ph.onons.E2l 
At half filling, the CDW state found at large V consists 
of doubly occupied sites without local moments, so that 
the spin order is strongly suppressed. The critical V c 
goes asymptotically to the value V c = 17/3 in both the 
weak- and strong-coupling limits, as predicted by mean- 
field theory. Between quarter and half filling, there is a 
mixture of singly and doubly occupied sites in the CDW 
phase, and the spin order is, in general, incommensurate. 
The total spin in the CDW phase is non-zero, i.e., charge 
order co-exists with ferromagnetism. As V is increased 
and the system goes into the CDW phase, the incom- 
mensurate spin correlations tend to be enhanced closer 
to quarter filling and suppressed closer to half filling. 
Note that this suppression of the spin moments within 
the CDW phase as the filling is increased is in contrast 
to the behavior in the homogeneous phase at V = and 
U ^$> t, for which the spin moment increases with increas- 
ing (n). 



An interesting subject for future study would be a more 
detailed comparison of the phase diag ram of the ladder 
system with that of the ID chain.li2lO For instance, the 
phase separation region in the ladder for large V/t is 
much smalleiEd than the \U\/t < 4 region found in ID. 
The origin of these differences in the charge-ordering pro- 
cess and the role of the non-zero ground-state spin remain 
to be clarified. 
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